Handedness of direct photons 
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The azimuthal asymmetry of direct photons originating from primary hard scatterings between partons is 
calculated. This can be accounted for by the inclusion of the color dipole orientation, which is sensitive to 
the rapid variation of the nuclear profile. To this end we introduce the dipole orientation within the saturation 
model of Golec-Biemat and Wiisthoff, while preserving all its features at the cross-section level. We show that 
the direct photon elliptic anisotropy v2 coming from this mechanism changes sign and becomes negative for 
peripheral collisions, albeit it is quite small for nuclear collisions at the RHIC energy. 



I. INTRODUCTION 

Direct photons can be a powerful probe of the underlying 
dynamics of the initial state of matter created in heavy ion col- 
lisions, since they interact with the medium only electromag- 
netically and therefore provide a baseline for the interpretation 
of jet-quenching models. There are several sources for direct 
photons, including prompt photons produced from initial hard 
scattering, thermal radiation from the hot medium and pho- 
tons induced by final state interactions with the medium. 

Unfortunately, the advantages of direct photons as a clean 
signature of the initial state of matter created in heavy ion col- 
lisions are offset by large backgrounds coming from hadronic 
decays, which should be extracted. The PHENIX collab- 
oration at RHIC has recently reported some results of the 
measurement of direct photon production [1, 2] , which has 
been also siibject of studies in several theoretical papers 

mwm. 

A novel mechanism which produces an azimuthal asym- 
metry coming from the reaction's initial conditions was in- 
troduced in Refs. IHSl- This is in contrast with the standard 
approaches where the azimuthal asymmetry is only associated 
with the properties of the medium created in the final state. In 
our approach, the main source of the azimuthal asymmetry 
originates from the sensitivity of parton multiple interactions 
to the steep variation of the nuclear density at the edge of the 
nuclei, which correlates with the color dipole orientation. In 
order to introduce a dependence on dipole orientation, we ex- 
tend the model of Golec-Biernat and Wiisthoff |9] for the total 
dipole cross section to the partial dipole-nucleon amplitude. 
To do that we assume that the two gluons in the Pomeron are 
not correlated. 



by a quark propagating and interacting with a target t (nu- 
cleon, t = N, or nucleus, f = A) at impact parameter b, as 
calculated from the diagrams in Fig. [1] can be written in the 
factorized form 1.3, 10.1 , 



d{lna)d^pd^b'"''^'"'' (27t) 

X (|)*^^ (a, n (a, r2 )/v an , ar2 , x) , 



(1) 



where p and a = py / p^ are the transverse and fractional 
light-cone (LC) momenta of the radiated photon and (^yg{a,r) 
is the LC distribution amplitude for the qy Fock component 
with transverse separation r. 
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Here Xij are the spinors of the initial and final quarks and 
Ko{x) is the modified Bessel function. The operators O have 
the form, 

d = inifa^ e* • (n X o) + ae* • (o x V) - i{2 - a) e* • V , (3) 

where e is the polarization vector of the photon, « is a unit 
vector along the projectile momentum, and V acts on rj- The 
parameter rug is the effective quark mass, which is in fact an 
infra-red cutoff parameter « 0.2 GeV. 

In equation ([T]i the effective partial amplitude 
Ft{b,ari,ar2,x) is a linear combination of qq dipole 
partial amplitudes at impact parameter b, 



II. PHOTON RADIATION IN THE COLOUR DIPOLE 
FORMALISM 

Radiation of direct photons in the target rest frame should 
be treated as electromagnetic bremsstrahlung by a quark in- 
teracting with the target. In the light-cone dipole approach the 
transverse momentum distribution of photon bremsstrahlung 



Ft(b,ari,ar2,x) = Im f Jb,ari,x) + f Jb,ar2,x 



fU{bMn-r2),x) 



(4) 



where x is Bjorken variable of the target gluons. The partial 
elastic amplitude f^^ can be written, in the eikonal form, in 
terms of the dipole elastic amplitude /^y of a qq dipole collid- 
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FIG. 1: Direct photons production in the target rest frame. The 
photon can be also radiated before the quark hits the target. Both 
diagrams are important. The dipole cross-section appears since the 
quark is displaced in impact parameter plane after radiation of pho- 
ton. The anti-quark enters after taking the complex conjugate of the 
amplitude, which is not shown here. 



ing with a proton at impact parameter b, 



1 



1-j I d^slmf^,{s,r,x)TA{b + s) 



l-cxp[- J f-slmf^^{s,r,x)TAib + s)]. 



(5) 

The hadronic cross section can be obtained by a convolution 
of the partonic cross section Eq. ([T]) with a proton structure 
function ^/(x,^) HI, 



da^pt yX) 
dxfd^prd^b 



da'^iqt yX) 
d{lna)d^PTd^D 



(6) 



where xp denotes the Feynman variable. We take the 
parametrization for the proton structure function given in 
Ref. d. 

We have recently shown that in this framework one can ob- 
tain a good description of the cross section for prompt pho- 
ton production data for prot on-proton (pp) collisions at RHIC 
and Tevatron energies flT]. Notice also that in contrast to 
the parton model, in this approach neither K-factor (NLO cor- 
rections), nor higher twist corrections are to be added. No 
quark-to-photon fragmentation function is needed either. In- 
deed, the phenomenological dipole cross section fitted to DIS 
data incorporates all perturbative and non-perturbative radia- 
tion contributions. Predictions for the LHC in the same frame- 
work are given in Ref. ifTsll . Comparison with the predictions 
of other approaches at the LHC can be found in Ref. iImIi . 



III. COLOUR DIPOLE ORIENTATION 

A colorless qq dipole is able to interact only due to the dif- 
ference between the impact parameters of q and q relative to 
the scattering center. If sis the impact parameter of the center 
of gravity of the dipole, and r is the transverse separation of 
the q and q, then the azimuthal angle of the radiated photons 
transverse momentum at a given impact parameter 5 correlates 
with the direction of s. In terms of the partial elastic amplitude 
fgg{s, r), it means that the vectors r and s are coiTelated. 
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FIG. 2: The relative distance of q and q from the center of gravity of 
qq dipole varies with the fractional light-cone momenta (3. 



One can see this in a simple example of a dipole interact- 
ing with a quark in Born approximation. The partial elastic 
amplitude reads, 

T f? r = ^ / d^qd^^i 

Jqis+r/2) _ ^iqis-r/l) 



Jq'iS+r/2) _iq'i^-7/2) 
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(7) 



where we introduced an effective gluon mass /j to take into 
account some nonperturbative effects. It is obvious from the 
above expression that the partial elastic dipole amplitude ex- 
poses a correlation between r and s, and the amplitude van- 
ishes when s - r = 0. In the above expression we assumed for 
the sake of simplicity that q and q have equal longitudinal 
momenta, i.e. they are equally distant from the dipole center 
of gravity. The general case of unequal sharing of the dipole 
momentum will be considered later. 

The Born amplitude is unrealistic, since it leads to an en- 
ergy independent dipole cross section agg{r,x). This dipole 
cross section has been well probed by measurements of the 
proton structure function at small Bjorken x at HERA, and 
was found to rise towards small x, with an x dependent steep- 
ness. 

The dipole elastic amplitude f^^ of a qq dipole colliding 
with a proton at impact parameter sis given by [|3D 

1 r d^qd^q' 
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where we defined a.; = \/ a.s{q^)(Xs{q'^) and J {x,q,q') is the 
generalized unintegrated gluon density (see below). The x de- 
pendence is implicit in the above expression. The fractional 
light-cone momenta of the quark and antiquark are denoted 
by p and 1 — p, respectively. It is obvious that the center of 
gravity of qq is closer to the fastest q or q, see Fig. (|2]l. The 
radiated photon takes away a fraction a of the quark momen- 
tum, see Fig. ([T]!. Therefore, for photon production, we have 



a 



(9) 
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FIG. 3: The partial elastic amplitude Im/^ (mb) of the qq dipole 
on a proton at impact parameter s as a function of dipole size r and 
angle 5 between .sand r for two values of a = 0, 1. We use a fixed 
value of X = 0.01 for all plots. 



Integrating over the vector s one can recover the dipole 



cross section a\:,-,(r,x 



4-K f (f-q 



{\-e-'^')aM')TM- (10) 



It is important to notice that the expression Eq. (|8]l also goes 
beyond the usual assumption that the dipole cross section 
is independent of the light-cone momentum sharing p. Al- 
though, the partial amplitude Eq. (O does depend on p, this 
dependence disappears after integration over impact parame- 
ter ? as shown in Eq. dTol i. 

The generalized unintegrated gluon density f {x,q,q') is 
related to the diagonal one by 



T{x,q,q' =q) = T {x,q). 



(11) 



The generalized unintegrated gluon density in Born approxi- 
mation takes the form. 



J {x,q,q') 
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FN{q-q')~F^^''\q,q' 



(12) 



where (A:) = (^Af|exp(/^-pi)|^Af) is the nucleon form fac- 
tor, and Fj^''\q,q') = (>PAr| exp[/^ • pi - iq' ■ ^2\\^n) is the 
so called two-quark nucleon form factor which can be cal- 
culated using the three valence quark nucleon wave function 
'Pa'(Pi,P2,P3)- 

For the dipole cross section we rely on the popular saturated 
shape |9] fitted to HERA data for F2{x^Q^). Assuming no 
correlation between the momenta q and q' inside the Pomeron 
aside from the Pomeron-proton form factor, we arrive at the 
following form of J {x,q,q') lH], 



3qo 

167i2a,. 
1 



X exp 



q'q''Rt>{x) 
Rl{x){q' + q'^) 



X exp[-4(^-^')V2], 



(13) 



where Oo = 23.03 mb, Roix) = 0.4 fm x (x/jco)"'^'*'* withxo = 
3.04 x 10^^ lUt]. We assume here that the Pomeron-proton 
form factor has the Gaussian form, F^{kj) = exp(— fc|./?^/2), 
so the slope of the pp elastic differential cross section is B'^j' = 
2Rjj + 2a'jpln{s/so), where ajp « 0.25 GeV"^ is the slope of 
the Pomeron trajectory, io = 1 GeV^. rJj sa (r^^^) /3 is the part 
of the slope of elastic cross section related to the Pomeron- 
proton form factor and (r^^J is the mean-square charge radius 
of the proton. 

Unfortunately, it is not possible to uniquely determine the 
unintegrated gluon density function from the available data. 
Nevertheless, the proposed form Eq. ST3[ seems to be a natural 
generalization which preserves the saturation properties of the 
diagonal part |3]. 

With this unintegrated gluon density the partial amplitude 
Eq. ([8]l can be calculated explicitly. 
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Rl 2Bel 



(14) 



where B,/(x)=/?2, + /;2(.x)/8. 

In Fig. (O we show the partial dipole ampUtude fqq{s,r) as 
a function of the dipole size r and the angle 5 between s and 
r, at various fixed values of s, for two values of a = 0, 1. One 



can see that for very small dipole sizes r the dipole orientation 
is not important. For very large dipole sizes r compared to the 
impact parameter s or very small values of s the dipole orienta- 
tion is also not present. It is important to note that the generic 
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FIG. 4: The anisotropy parameter {b,p,a) as function of 
calculated at a = 1 for different impact parameters b and energies: 
^ = 200GeV (solid, b = 0.2, 0.4, 0.6, 1 fm), and Jl = 5500GeV 
(dashed, b = 0.6, 1 fm). The plot is taken from Ref. [3j]. 



feature of the partial dipole ampHtude, e. g. its maximum and 
minimum pattern, changes with a. However, it is not obvious 
a priori how the convolution with the proton structure function 
Eq. (|6]l, which leads to a sum over all different configurations 
of a and on top of that the convolution between the partial 
dipole amplitude and the nuclear profile Eq. (|5]l, which leads 
to a even more complicated angle mixing, gives rise to a final 
azimuthal asymmetry. 



IV. AZIMUTHAL ASYMMETRY 

The main source of azimuthal asymmetry in the amplitude 
(|5]l is the interplay between multiple rescattering and the shape 
of the physical system. The key function which describes the 
effect of multiple interactions is the eikonal exponential in 
Eq. (|5]l, while the information about the shape of the system 
is incorporated through a convolution of the impact parameter 
dependent partial elastic amplitude and the nuclear thickness 
function. Notice that the initial space-time asymmetry gets 
translated into a momentum space anisotropy by the double 
Fourier transform in Eq. ([T]i. 

The azimuthal asymmetry of prompt photon production, re- 
sulting from parton-nucleus (qt) or proton-nucleus (pt) colli- 
sions for t ~ N,A, is defined as the second order Fourier co- 
efficients in a Fourier expansion of the azimuthal dependence 
of a single-particle spectra Eq. ([T]i around the beam direction. 



{PT,b,a) 



where the angle (|) is defined with respect to the reaction plane. 
In the same fashion, the azimuthal asymmetry of photon yield 
from collisions of two nucleus Ai and A2 at impact parameter 
B is defined as 



/!^„Qf^COs(2^) gN _ 

.^gcos(20O^^^$r^7,,(/;.) 

Art (nr^ .d^KpM-^jX) 



dxfd^pjd^bi 



d b 1 t^Ta, (/JiJ 



dxfd^PTd^b2 



(16) 



(15) 



where we used the notation b2 = b + B, bi = b (b is the im- 
pact parameter of the pAi collision) and the angle ©i (©2) is 
the angle between the vectors bi( ^2) and B, respectively. The 
medium modification of nucleon structure functions in our in- 
terested range of pr is less than 20% and is ignored in the 
above expression. 

The only external input in our approach is the nuclear pro- 
file. First, we take a popular Woods-Saxon (WS) profile, with 
a nuclear radius Ra —6.5 fm and a surface thickness ^ = 0.54 
fm, for PbH-Pb colhsions lUsll . 

In Fig.m we show examples of azimuthal anisotropy from 
quark-nucleon collisions radiating a photon vf^lb, p, a), with 
a = 1 and at different impact parameters and energies. The 
results show that the anisotropy of the dipole interaction rises 
with impact parameter, reaching rather large values. As func- 
tion of the transverse momentum of the radiated photons, 
vf^ {b,p,a) vanishes at large pj. Such a behavior could be 
anticipated, since the interaction of vanishingly small dipoles 
responsible for large p is not sensitive to the dipole orienta- 
tion. 

In Fig. (|5]), we show the calculated values of defined 
in Eq. ( fTSl l. for fixed a = 1, at various q(p)A collision impact 
parameters b for the RHIC energy \/s — 200 GeV at midra- 
pidities. If the nuclear profile function was constant, then the 
convolution between the nuclear profile and the dipole orien- 
tation, defined in Eq. (|5]l, would be trivial, and v'2' becomes 
then identically zero. Therefore, the main source of azimuthal 
anisotropy is not present for central collisions where the cor- 
relation between nuclear profile and dipole orientation is min- 
imal. This can be seen in Fig. (|5]l, where a pronounced elliptic 
anisotropy is observed for collisions with impact parameters 
close to the nuclear radius Ra, where the nuclear profile under- 
goes rapid changes. Therefore, the important parameter which 
controls the elliptic asymmetry in this mechanism is \b — RA\ 
13]. 

It is important to notice that V2 is suppressed an order of 

magnitude compared to Vj^. At first glance this might look 
strange, since the quark interacts with nucleons anyway. How- 
ever, a quark propagating through a nucleus interacts with dif- 
ferent nucleons located at different azimuthal angles relative 
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cellation would be exact if the nuclear profile function Ta {b) 

aA 

were constant. We have a nonzero, but small Vj only due 
to the variation of Ta with b. Going from q{p)A to AA colli- 
sions, see Fig. (|5]l, the elliptic asymmetry is further reduced. 
The main reason is that the integrand in Eq. ( fT6] l gets con- 
tributions only from semi-peripheral pA collisions where our 
mechanism is at work, and most of the integral over b does 
not contribute. This significantly dilutes the signal. 



V. ON THE SIGN OF V2 

In Fig. (O, we showed that the general behaviour of dipole 
amplitude orientation, e.g. its maximum and minimum pat- 
tern changes with the parameter p which defines the relative 
position of the center of gravity of the dipole from the q and q, 
see Fig. (|2]i. Here, we take a heuristic approach, and explore 
a possible link between the peculiar behaviour of V2 shown in 
Fig. (|5]l with the dipole orientation introduced in Eq. (O. 

Let us assume for sake of argument that the parameter p is 
independent of a and assume that q and q have equal longitu- 
dinal momenta namely p = 1/2. This corresponds to a partic- 
ular configuration in which the dipole amplitude is symmetric 
under ? — r, see Fig. (O. In principle, this configuration is 
kinematically less probable for direct photon production (in 
contrast to DIS) since it corresponds to a = via Eq. (|9]l. No- 
tice that although we take a fixed p = 1/2 in the dipole am- 
pUtude Eq. (O, the LC distribution of the projectile quark yq 
fluctuation Eq. ^ still depends on a and also the transverse 
dipole size is ar and varies with a, see Eq. ([T). 

We repeat the computation of the azimuthal asymmetry V2 
of prompt photons in the same way as discussed in the previ- 
ous section. For example, in Fig. (|6]l, we show the anisotropy 



asymmetry 



and 



qA 



V2 at a = 1 for various impact parame- 
ters as in Figs. (1415b . Comparing with the results presented in 
the previous section, it is seen that although the order of mag- 
nitude of V2 is the same in both cases, now the sign of V2 does 
not change at higher pj and remains positive. This indicates 
that the sign of V2 in this mechanism is related to the dipole 
orientation via the parameter p. 

Notice also that the sign behaviour of the prompt photon V2 
for AA collisions at higher pj is also present for both qN and 
qA collisions. 



VI. CONCLUDING REMARKS 



FIG. 5: The impact parameter dependence of prompt photon az- 
imuthal asymmetry, for q-l-Pb (at a = 1), p+Pb and Pb+Pb collisions 
for RHIC energy at midrapidities, with the Woods-Saxon (WS) nu- 
clear profile. As example, we have also shown in the first panel, the 
prompt photon elliptic anisotropy from q+pb collisions for the hard 
sphere (HS) nuclear profile. The plot is taken from Ref. 0]. 



to the quark trajectory. Their contributions to tend to can- 
cel each other, restoring the azimuthal symmetry. Such can- 



The azimuthal elliptic asymmetry V2 observed in heavy ion 
collisions is usually associated with properties of the medium 
created in the final state. We introduced a novel mechanism 
which relates this azimuthal asymmetry to the colour dipole 
orientation. To this end, we proposed a model generalizing 
the unintegrated gluon density fitted to data for the proton 
structure function to an off-diagonal unintegrated gluon dis- 
tribution. 

We showed that the azimuthal asymmetry V2 of prompt pho- 
tons changes sign and becomes negative for peripheral colli- 
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FIG. 6: Upper panel: The anisotropy parameter I'j {b, p, a) as func- 
tion of pT for different impact parameters b for RHIC and LHC ener- 
gies. Down panel: The anisotropy parameter v^{b,p,a) as function 
of px for various impact parameter b for RHIC energy at midrapid- 
ity. Similar to Fig. {5) we show the results for both the WS and the 
HS nuclear profiles. In both plots the results calculated at a fixed 
P = 1/2 in the dipole amplitude and a = 1. 



sions. Although this behaviour seems to be robust for the con- 
sidered range of pr, there is some uncertainty on the mag- 
nitude of the v2 coming from this mechanism. The shape 
of the tail of nuclear profile is important in this mechanism, 
since it significantly affects the results. To highlight this point, 
in Figs. ( I5I6I ). we have also shown the azimuthal asymmetry 
from quark-nucleus collisions for the hard sphere (HS) 
nuclear profile. By comparing with the results from the WS 
nuclear profile for the same setting, one may conclude that the 
maximum uncertainty in this mechanism can be as big as an 
order of magnitude. Unfortunately the tail of all available nu- 
clear profile parametrizations is less reliable and obtained by 
a simple extrapolation |15]. This is also due to the fact that 
the neutron distribution, which may be more important on the 
periphery, cannot be properly accounted for by electron scat- 
tering data. Another source of uncertainty in this approach is 
due to the fact that the off-diagonal part of the unintegrated 
gluon density cannot be uniquely defined from the current ex- 
perimental data. 

In order to see if this mechanism is relevant to heavy ion 
collisions, it would be of great interest to calculate the az- 
imuthal asymmetry for gluon radiation and hadron produc- 
tion at RHIC. This mechanism might also contribute to the 
azimuthal asymmetry in DIS and in the production of dilep- 
tons. We plan to report on some of these problems in the near 
future. 



Acknowledgments 



One of us (AHR) would like to thank the organizers of "// 
Latin American Workshop on High Energy Phenomenology" 
for the kind invitation and for the very interesting and stim- 
ulating meeting. This work was supported in part by Coni- 
cyt (Chile) Programa Bicentenario PSD-9 1-2006, by Fonde- 
cyt (Chile) grants 1070517 and 1050589, and by DFG (Ger- 
many) grant PI182/3-1. 



[1] PHENIX Collaboration, arXiv:0705.17Il 

[2] PHENIX Collaboration, Phys. Rev .Lett. 96, 032302 (2006). 

[3] B. Z. Kopeliovich, H. J. Pirner, A. H. Rezaeian, Ivan Schmidt, 

Phys. Rev. D77, 034011 (2008)|arXiv:0711.3010|. 
[4] B. Z. Kopeliovich, A. H. Rezaeian, Ivan Schmidt, 

ar Xiv:07 12.2829 Accepted for publication in Nucl. Phys. 

A. 

[5] S. Turbide, C. Gale and R. J. Fries, Phys. Rev. Lett. 96, 032303 
(2006). 

[6] S. Turbide, C. Gale, E. Frodermann, U. Heinz, Phys. Rev. C77, 

024909 (2008) |arXiv:0712.0732|. 
[7] R. Chatterjee, E. S. Frodermann, U. W. Heinz, D. K. Srivastava, 

Phys .Rev. Lett. 96, 202302 (2006)^^ 

[8] F M. Liu and K. Werner, arXiv: 0712.3619l 



[9] K. Golec-Biernat and M. Wusthoff, Phys. Rev. D59, 014017 
(1999). 

[10] B.Z. Kopeliovich, A. Schaefer and A.V. Tarasov, Phys. Rev. 

€59 (1999) 1609 
[11] B. Z. Kopeliovich, A. H. Rezaeian, H. J. Pirner and Ivan 

Schmidt, Phys. Lett. B653, 210 (2007). 
[12] SMC Collaboration, Phys. Rev. D58, 112001 (1998). 
[13] A. H. Rezaeian, B. Z. Kopeliovich, H. J. Pirner and Ivan 

Schmidt, arXiv:0707.2040 
[14] S. Abreu et al, arXiv: 071 1.0974 

[15] A. Bohr and B. R. Mottelson, Nuclear Structure (Benjamin, 
New York, 1969); C. W. De Jager, H. De Vries, C. De Vries, 
Atom. Data Nucl. Data Tables 36, 496 (1987). 



